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Abstract. We study the effect of quenched disorder on the zero-range process (ZRP), a 
system of interacting particles undergoing biased hopping on a one-dimensional periodic 
lattice, with the disorder entering through random capacities of sites. In the usual ZRP, 
sites can accommodate an arbitrary number of particles, and for a class of hopping rates and 
high enough density, the steady state exhibits a condensate which holds a finite fraction of 
the total number of particles. The sites of the disordered zero-range process considered here 
have finite capacities chosen randomly from the Pareto distribution. From the exact steady 
state measure of the model, we identify the conditions for condensate formation, in terms of 
parameters that involve both interactions (through the hop rates) and randomness (through 
the distribution of the site capacities). Our predictions are supported by results obtained 
from a direct numerical sampling of the steady state and from Monte Carlo simulations. 
Our study reveals that for a given realization of disorder, the condensate can relocate on 
the subset of sites with largest capacities. We also study sample-to-sample variation of 
the critical density required to observe condensation, and show that the corresponding 
distribution obeys scaling, and has a Gaussian or a Levy-stable form depending on the 
values of the relevant parameters. 
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1. Introduction and model 


Quenched disorder can strongly affect both static and time-dependent properties of statistical 
systems. Of particular interest is the case of driven systems in which a dynamics that violates 
detailed balance leads the system to a nonequilibrium stationary state (NESS) that cannot 
be described within the purview of the Boltzmann-Gibbs equilibrium statistical mechanics. 
In this work, we explore effects of quenched disorder by analyzing a one-dimensional model 
of a disordered nonequilibrium system, namely, a disordered zero-range process. The model 
is a modification of the well-studied zero-range process (ZRP), a lattice model of interacting 
particles evolving in presence of an external drive, with no limit on the capacity of each site 
to hold any number of particles CHI. Here, we study a disordered model introduced in [5][6], 
in which the capacity of each site is a randomly chosen finite number. Knowledge of the 
exact steady state of this model allows us to unveil and understand the physical effects that 
result from quenched disorder. We are primarily interested in the possible occurrence of a 
condensate, which is a quintessential feature of the ZRP, as discussed below. 

On a one-dimensional periodic lattice, the ZRP dynamics involves particles undergoing 
stochastic hopping between the lattice sites mat For a system of L sites and N 
indistinguishable unit-mass particles, a unit time step of the dynamics comprises L sequential 
moves, in each of which a particle hops out of a random site i; * = 1, 2,..., L, with occupancy 
rip rii > 0, with a specified hop rate Ui(rii), and moves to site i + 1. The particle density 
is p = N/L. The forward-biased hopping of a particle from a site to only its right neighbor 
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incorporates the effect of an external driving field on the particles. Evidently, the dynamics 
conserves the total number of particles in the system. While there is no interaction between 
particles on different sites, that between particles on the same site may be modeled through 
the dependence of the hop rate of a particle from the site on its occupancy. Remarkably, 
the NESS measure of configurations in the ZRP can be found exactly for any choice of hop 
rates and in any spatial dimension mat The homogeneous ZRP is defined by having the 
same functional form of the hop rate for all sites. 

The phenomenon of real-space condensation that can occur in the steady state of the 
ZRP involves a finite fraction of particles accumulating on a single site, thereby forming a 
macroscopic condensate whose mass increases with increasing density. In the case of the 
homogeneous ZRP, a possible choice of the hop rate is u(ri) = 1 + b/n, where b > 0 is a 
finite constant. Such a form of the hop rate implies an effective attraction between particles 
on the same site. For such a choice, it is known that for b > 2, the model in the NESS 
exhibits a transition to a condensate phase at a critical value of the particle density given 
by p c = 1/(6 — 2) [3], For densities p < p c , the system is in a fluid phase characterized by 
an occupancy of order unity on every site and a single-site occupancy distribution p{n) that 
decays exponentially for large n. At the transition point p = p c , the distribution decays 
asymptotically as a power-law, corresponding to a critical fluid. Above p c , the critical fluid 
coexists with a macroscopic aggregate (the “condensate”), so that in addition to a power- 
law part, the distribution p{n) has a sharp peak around n = (p — p c )L that represents the 
condensate. The ZRP has been invoked to model condensation in a number of contexts, 
e.g., clustering of particles in shaken granular systems [7], jams in traffic flows (8], wealth 
condensation in macroeconomics [9j, and other systems. 

Driven diffusive systems constitute a class of stochastically evolving interacting particle 
systems typified by a spreading of density fluctuations with a systematic drift in addition to 
a diffusive motion [T0I - H2] . At long times, these systems relax to a NESS in which a steady 
current of particles flows through the system. The forward-biased ZRP described above is 
an example of a driven diffusive system, and can also be mapped to another paradigmatic 
and extensively studied model in this class, namely, the asymmetric simple exclusion process 
(ASEP). On a one-dimensional periodic lattice, the ASEP involves indistinguishable hard¬ 
core particles undergoing biased hopping to empty nearest-neighbor sites. The mapping 
between the ZRP and the ASEP consists in interpreting sites (respectively, particles) in the 
former as particles (respectively, empty sites) in the latter [3]. 

Quenched disorder in driven diffusive systems has been studied over the years in 
several types of systems in this class 0. Both particle-wise disorder, in which different 
particles have different time-independent hop rates P3HI8], and space-wise disorder, with 
time-independent hop rates that are randomly distributed in space [5ll6L fT9l - [24] . have been 
considered. Note that the mapping between the ASEP and the ZRP mentioned above 
transforms particle-wise disorder in the ASEP into space-wise disorder in the corresponding 
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ZRP (THUS]. Other studies of quenched disorder in driven systems include a disordered 
ASEP with particle non-conservation, in which randomly chosen sites do not conserve particle 
number [25], a ZRP on inhomogeneous networks, in which particles hop between nodes of 
a network with one node of degree much higher than a typical degree [26], a ZRP with 
quenched disorder in the particle interaction, implemented through a small perturbation 
of a generic class of hop rates [27, 28], and a ZRP involving an interplay between on-site 
interaction and diffusion disorder [29]. 

In contrast to the above mentioned studies of quenched disorder, in the model under 
consideration here, disorder is assigned to the capacities of sites; the capacity is a random 
variable that restricts the number of particles a site can accommodate. This model was 
introduced in [HIE] as the disordered drop-push process (DDPP) to study transport of 
carriers trapped in local regions of space, and is a generalization of the uniform drop-push 
process [30][3l]- As pointed out in [2], the drop-push process is actually a special case of the 
ZRP, with infinite hopping rates out of sites in which the occupancy exceeds the capacity. 
We thus prefer to refer to the model as the random capacity zero-range process (RC-ZRP). 
The fact that capacities are finite and random has a strong and an essential influence on the 
ZRP steady-state dynamics, as we discuss below. 

In this paper, we consider capacities chosen independently for every site from a common 
distribution with power-law tails. For a one-dimensional periodic lattice with L sites, every 
site i has a capacity U chosen independently from the Pareto distribution: 



Note that for a given realization {h} of the disorder, the system can accommodate at most 

L 

N = V /■ 

1 v max — / v 



A new feature, namely, a dynamical cascade effect, emerges owing to sites having 
restricted capacities in the RC-ZRP [5] El EDI El]. Consider a particle hopping out of a 
random site i that has occupancy 0 < n* < Z* with hop rate tq(nj|Zj), and moving to site 
i + 1. If the site i + 1 is already full, a particle from this site gets pushed further right, and 
so on, leading to a sequence of adjacent-site hops that continues until a particle hops into a 
site ( i + m , say) that was not fully occupied earlier (i.e., n i+m < li +m ). Note that a unit time 
step corresponds to L updates at randomly chosen sites, where each update may involve 
several particle hops out of fully filled sites. Thus, restricted capacities lead to a cascade of 
particle transfers through filled sites, explaining the nomenclature “drop-push process” used 
in [30][3T]. Figure Q] shows a schematic view of the RC-ZRP. 

In this paper, we ask: Does restricting the ZRP site capacity, as in the RC-ZRP, still 
allow for the formation of a condensate, and if so, under what conditions? Let us consider a 

| In simulations of the RC-ZRP reported later in the paper, the capacity of a site is taken as the largest 
integer not exceeding a real number drawn from the distribution ©• 
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Figure 1. Schematic diagram of the random capacity zero-range process, showing a typical 
configuration and possible particle hops. 


form for the hop rate Ui{n\li) that promotes condensate formation in the homogeneous ZRP, 
namely, Ui(n > 0) = 1 + b/n V i, with b > 0. The hop rate is taken to have this form for 
occupancies in the range 0 < n < li, and to be infinite for occupancies larger than Zj |2j, 
ensuring an immediate movement of a particle from a filled site to one that is not full: 

{ 0 for n = 0, 

1 + y for 0 <n <k, (2) 

oo for n > k. 

Note that in contrast to earlier studies of the model in M which considered general hop 
rates, the choice ([2]) has the possibility of supporting condensate formation. In this paper, we 
demonstrate on the basis of exact analytical and simulation results that an interplay of the 
capacity distribution ((Tj) with the hop rate ([2]) can indeed lead to condensate formation, and 
derive the conditions for this to happen. It should be noted that the steady state equal-time 
properties reported in this paper hold not just for the considered case of biased hopping of 
particles from a site to its right neighbor, but also for unbiased hopping to the left and to the 
right neighboring site. This is because, as we discuss in Section [2l the steady state measure 
of the RC-ZRP is the same in the two cases. However, unequal-time properties in the steady 
state, for instance, relocation dynamics of the condensate, will be different for biased and 
unbiased hopping. 

The ZRP in which sites have bounded capacities was addressed recently in [32]. Unlike 
our model, the capacity was taken to be the same for all sites, and quenched disorder was 
introduced through site-dependent and particle-dependent hop-rates, leading to dynamical 
blocking that causes slow relaxation to steady state. 
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The layout of the paper is as follows. In Section [2l we discuss the steady state of the 
RC-ZRP, based on which we derive in Section [3] the conditions to obtain condensation in 
the model. In Section |4j we confirm our predictions by a direct numerical sampling of the 
steady state and by performing Monte Carlo simulations of the steady-state dynamics. We 
also discuss the relocation dynamics of the condensate. The paper ends with conclusions 
in Section [5] while the Appendix summarizes some relevant features associated with the 
capacity distribution. 


2. Stationary state 


The RC-ZRP relaxes at long times to a current-carrying nonequilibrium stationary state. 
Using the condition of pairwise balance ra , as in OT, the steady state measure of 
configurations may be found. For a given realization {l A of the disorder and a given total 
number of particles N < A r max , it has a factorized form | 


'P({n i }\{l i },N) oc Wfi{ni\li)5 ( ) , 


(3) 


2=1 


, 2 — 1 


where 5(m,n) is the Kronecker Delta function, while the single-site factors /,.(n|Zj) equal 
unity for n = 0 and are given for n > 0 by 

-l 


fi{n\k) = 


J Ui(m\li) 


771=1 


I Ufe+n+i) for 0 < n < i 

_ ) r(n+i)r(6+i) iUI u ^ u — 

0 for n > k , 


(4) 

(5) 


where r(x) is the Gamma function. Here, in arriving at the second equation, we have used 

Eq. ©. 

Equation (J3]) is the measure of configurations within the canonical ensemble. In the 
thermodynamic limit N —> oo ,L —>- oo, keeping the overall particle density p = N/L fixed, 
we use an equivalent grand canonical ensemble description of the steady state. In such a 
description, the total number of particles is allowed to fluctuate, and a fugacity z fixes the 
average number of particles to equal N. The steady state measure of configurations within 
the grand canonical ensemble is given by 

L 

Prob({ni}|{Zj}) oc '\\pi{ni\k), (6) 

2=1 


§ The measure (j3j) holds in any spatial dimension, for any choice of the hop rate, and for any rules of particle 
transfer, either biased or unbiased, between sites. In case of unbiased transfer, a case not addressed here, 
detailed balance holds, and the steady state is an equilibrium state. 
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where Pi{n\lj) is the single-site occupancy distribution, namely, the probability for the ith 
site to have 0 < n < U particles: 


Pi(n\li) 


z n fi(n\li ) 

Fi{z\k) ’ 


(7) 


with Fi(z\U) ensuring normalization of Pi(n\li): 


F i {z\l i ) = l + '£z n f i {n\l i 


( 8 ) 


71=1 


Here, the fugacity z satisfies Y2i=i™i = N, where nj = dinFi(z\li)/dinz is the average 
occupancy at the ith site, with the average taken with respect to the single-site probability 
dH) . We finally get 




zF'Mk) 


= P . 


(9) 


Lti Fi(Ak) 

where prime denotes differentiation with respect to z. 

As discussed in the Appendix, the maximum number of particles N max that can be 
accommodated in the system scales with L as N mgbX ~ L 1//q for a < 1, and as N max = aL for 
a > 1, where a = f dl IP{l ) is finite. It then follows that in the latter case, the density p 
in Eq. (J9J) has a maximum allowed finite value equal to a, while for a < 1, p diverging with 
L as L l/a- 

1 can be arbitrarily large. For a given value of z, evaluating numerically the left 
hand side of Eq. (j9j) for a given realization of disorder, and then averaging with respect to 
disorder, we show in Fig. [2]the disorder-averaged p, denoted by (p), as a function of z for 
three different system sizes L[jJ. 


|| Here and in the rest of the paper, angular brackets will be used to denote averaging with respect to disorder 
realizations. 
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Figure 2. Disorder-averaged density (p) versus fugacity z, computed using Eq. |9]) for 
b = 2.5 and a = 0.5. The system sizes are marked in the figure. Here, Ay, is the number of 
disorder realizations over which the data have been averaged. 


The occupancy distribution for the full system is defined as the probability that a 
randomly chosen site has n particles. This is possible only if the site capacity is equal to or 
larger than n, so that the distribution has the form 


/»oo 

pin) oc z n f{n ) / dl P(l). 

(10) 

For large n 1, Eq. ((5]) combined with the approximation F{b + n 

n b + Oin^ 1 ) yields the following asymptotic behavior: 

+ l)/r(n + 1) « 

exp (—n/n*) 

Pi{n\li) x h 

n° 

(11) 

. exp(— n/n*) 

K n b + a > 

(12) 

where the characteristic occupancy n * is given by 


n * = —1/ ln(z). 

(13) 


The RC-ZRP in the NESS supports a steady current of particles through the system. 
Within the grand canonical ensemble, an exact expression for the average steady-state 
current for a given realization of disorder was derived in [6]. We briefly summarize the 
derivation here. From the dynamics, it is evident that all hops contributing to the current 
across the bond (* — 1 ,i) for which site i is completely full also contribute to so 

that one has the recursion 

k 

Ji,i+1 Pi (^i | ^i) Ji—l,i T ^ ^ \li)Pi {j^i |^i) j 

Ui = 1 

8 


L=4,N h =500 + 

L=8,N h =250 
L=16,N h =250 + 










4 < 


(14) 





where the second term on the right hand side is due to hops originating from site i. Now, 
steady state implies that all the bond currents are equal, J iji+ 1 = = ... = J 0 , so that 

the above equation yields 


j Y2n i =l Ui ( ni \l i )P i ( n i\h) 

J ° = -1-777m- = 

1 Pi\J>i \h) 


(15) 


where we have used Eqs. (J4j) and ([7]) . The steady-state current depends through the fugacity 
z on the overall particle density p and the number of sites of different capacities in a given 
realization of disorder (see Eq. ()9|) ), thus becoming a function of b and a. Note that the 
expression (jT5l) for the average steady-state current in terms of the fugacity 2 is the same as 
for the usual ZRP [3]. 

In the following section, we address the issue of condensate formation in the RC-ZRP. 


3. Condensate formation 

In this section, we turn to the conditions for condensate formation in the RC-ZRP. We also 
study the distribution of the critical density p c to observe condensation, and its scaling as a 
function of L. We begin by summarizing the main questions and the results before getting 
to the details of the derivation. 

It is useful to first recall the known scenario of condensation in the customary 
homogeneous ZRP [3]. Since there is no restriction on site capacities, there is no difficulty in 
accommodating O(L) particles on any one site. An essential requirement for condensation is 
the existence of a hnite critical value of the average site occupancy in the limit the fugacity 
z attains its maximum possible value £ max . Below this critical value of the density, all sites 
have the same average occupancy equal to the overall density of particles in the system. 
Above the critical value, the average occupancy of all but one site has the critical value; the 
excess particles that form a hnite fraction of the total number of particles are accommodated 
on a single randomly-chosen site. 

In this backdrop, it is a priori not apparent whether and when such a scenario of 
condensation holds in the RC-ZRP in which sites have restricted capacities. To address the 
issue, we argue as follows. 

(i) A necessary condition for condensate formation is that at least one site be able to 
accommodate O(L) particles. In view of sites having restricted capacities in the RC- 
ZRP, the candidate for a site that can accommodate a macroscopic number of particles is 
the one with the largest capacity. Then, if the largest capacity / max = Max[Zi, / 2 ,..., II] 
has the scaling / max ~ L e , we need 6 to be larger than unity for condensate formation. 

(ii) Additionally, we require that the average site occupancy has a hnite value, denoted by 

p ci as * ^ ^ maX ' 
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(iii) When conditions (i) and (ii) are fulfilled, condensate formation is possible at high enough 
density p > p c . The critical density p c depends on the realization of disorder, and the 
question arises as to how the form of the disorder-induced distribution Prob(p c ) of p c 
depends on the relevant parameters. 

Whether conditions (i) and (ii) above would hold depends on parameters that 
characterize the probability distribution of capacities (Eq. (JTj)) and the form of the hop 
rate (Eq. ((2])), namely, the exponents a and b. The results are as follows: 

(i) The exponent 6 is given by 1/a, so the site with the largest capacity can accommodate 
O(L) particles provided that a < 1. 

(ii) The average site occupancy remains finite as z —> £ max so long as one has b + a > 2. 
Combining the last two points, we thus arrive at the following conditions for condensate 
formation in the RC-ZRP: 

b + a > 2, and a < 1 (Conditions to obtain condensation). (16) 

(iii) The distribution Prob(p c ) is a Gaussian for b > (4 — a)/ 2, while it is a Levy-stable 
distribution for 1 < b < (4 — a)/2. 

These results are summarized in Fig. EJ which shows the regime for condensate formation 
in the a — b plane, and also the forms of Prob(p c ) in different regions. 

The issue of sample-dependence of phase transitions was studied numerically in [21] for 
the one-dimensional ASEP with open boundaries and quenched-disordered hopping rates. 
In the present work, we are able to determine the analytic forms for the distribution of the 
critical density because of the product form of the steady state measure in the RC-ZRP with 
periodic boundary conditions. 

In the limit a —> 0, the capacities become infinitely large, and the RC-ZRP dynamics 
becomes similar to the dynamics of the homogeneous ZRP. In this limit, the condition to 
observe condensation becomes the requirement b > 2, a result known for the homogeneous 
ZRP [3]. 

We now proceed to a derivation and a more detailed discussion of our results. 
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Figure 3. In the a — b plane, the dashed regions are those for which the system exhibits a 
condensate; the corresponding forms of the p c distribution are also indicated. 


3.1. The largest capacity Z max 

In order to accommodate the condensate, it is necessary that Z max grows sufficiently rapidly 
with L , namely, as L 9 , with 6 > 1. For the Pareto distribution of Eq. ([I]), it is known that 
Zmax scales as L 1 /", see the Appendix. Thus, condensate formation is possible provided that 
a < 1 . 

As discussed in the Appendix, not just the site with capacity Z max , but in fact several 
sites have capacities of order L 1 /". This feature allows the condensate to form and relocate 
in time on other sites that belong to this subset. The numerical studies reported in Section 
[4] bear this out. 

3.2. Average site occupancy 

Let us now investigate the behavior of the average site occupancy as z — >• z max , in the regime 
a < 1 as required to accommodate a putative condensate (Section I3.1H . The quantity z max 

l max 

is obtained by requiring the convergence of the series Fi 1 (z) = 1 + Y2 z n fi 1 (n) in the limit 

n =1 


li 













Figure 4. G(l) as a function of l for three representative values of b , obtained by evaluating 
numerically Eq. m- The black lines stand for the asymptotic behavior, Eq. m- 


L —* oo, where i\ is the label for the site with the largest capacity. We obtain z msix = 1 as 
the radius of convergence of the series. The average site occupancy being a monotonically 
increasing function of z has a maximum allowed value in the limit z —» z max , given by 

/({;<})= lim ljr zFmi) 


Z —>-Zmax L ' ^ 


7 £g(0; 


G ^ = 


i= 1 


/ , (17) 

i + Eti fMi) 

which defines a characteristic density p * for every disorder realization. Evidently, the density 
p*({h}) is given by a sum of L i.i.d. random variables G{li ); i = 1, 2,..., L. A finite value 
of p* implies condensate formation for p > p*, so that p* coincides with the critical density 
p c to obtain condensation. From Eq. (117)) . we get the corresponding disorder-averaged 
characteristic density as 


/ OO 

dl P{l)G{l). (18) 

The behavior of p * is governed by two parameters, namely, the exponent a that 
characterizes the probability distribution of the capacities, and the exponent b that 
characterizes the hop rate. Let us ask for the condition on the allowed range of values 
of b and a that leads to a finite (p*). Using Eq. ([5]), the function G(l ) in Eq. (1171) can be 
expressed as 


G{1) 


1 _ (6 - !)(/ + 1 ) _ 

h ~ 2 {b - 2) (l - T(6 + l + 1)/[T(6)T(Z + 2)]) ’ 


(19) 
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implying that G{1) diverges for particular values 6 = 1,2. In the asymptotic regime l 1, 
using r(6 + / + l)/r(Z + 2) l b ~ l + 0(l b ~ 2 ) } we get 


G(l > 1) « 


(1 ~b)l 
2—b 

( b-l)T(b)l 2 ~ b 
2—b 
1 

b—2 


for 0 < 6 < 1, 
for 1 < 6 < 2, 
for 6 > 2. 


( 20 ) 


Figure |4] illustrates that the asymptotic form, Eq. (]20]) . is a good approximation to the exact 
expression, Eq. (TT9|) . Now, Eq. (TT8|) gives 

/»oo 

(p*)«r(Z*)+/ dl P(l)G(l), (21) 

where l* is chosen such that for / > l*, the function G{1) is well approximated by its 
asymptotic behavior, Eq. ([20)1 . for large l. The value of l * depends on 6; for example, 
from Fig. [I] one may choose l* = 100 for 6 = 0.5 and l* = 1000 for 6 = 1.5. In Eq. (Ell) , the 
finite constant r(l*) is the value of the integral fl dl P(l)G(l). By analyzing the integral in 
Eq. (El]) , one then concludes that requiring (p*) to be finite leads to the following conditions: 


a > 1 for 0 < 6 < 1, 

6 + a > 2 for 1 < 6 < 2, (22) 

a > 0 for 6 > 2. 


The above conditions may be combined into the single condition 

6 + a > 2 for (p*) to be finite. (23) 


This is to be contrasted with the condition 6 > 2 in the homogeneous ZRP for the average 
site occupancy to be finite as z —* £ max - 

At the critical point, when z —> z max , the characteristic occupancy n* in Eq. (Il3l) 
diverges, and we find from Eqs. ([IT]) and (TT2]) the occupancy distributions for n 1 to obey 


lirn Pi(n\li) oc —r, 

Z —^max Tl° 


lim p(n) oc 

Z->2max 


1 

n b+a ' 


(24) 

(25) 


3.3. Distribution of the characteristic density p * 

For given values of L, ex, and 6, one may obtain the density p*({h}) for different disorder 
realizations {U} by using Eq. ([TTj) . Let us denote the corresponding distribution as 
Prob(p*). One may deduce the form of Prob(p*) by invoking the well-known theory of stable 
distributions, which concerns the sum S L = °f a number L of mutually independent 

random variables {Aj; i = 1,2, ...,L} having a common distribution [33ff35] . When this 
common distribution has a power-law tail decaying as |A01 then, for 0 < cx < 2, the 
limiting distribution for Sl as L —» oo converges in form to a stable distribution that has 
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Figure 5. Illustrating the validity of the behavior (l26l) of the probability distribution 
Prob(G) for three representative values of b at a = 0.375. The values of G are computed 
using Eq. m- The data for b = 1.5 have been scaled up by a factor of 300 for convenience 
of display. The black lines denote analytical predictions, namely, (i) for b < 2, a power-law 
decay with exponent (1 + a) for 0 < b < 1, and exponent (1 + v) for 1 < b < 2, and (ii) for 
b > 2, a unit step function. 


a tail decaying as |<S'l| _ 1_Q! - For a > 2, on the other hand, the distribution converges to a 
Gaussian distribution. 

From the large-/ behavior of G(l) in Eq. (12011 . we deduce the following tail behavior of 
Prob(G): 

{ G~ 1_ “ for 0 < b < 1, 

G~ l ~ v ] v = aj (2 - b) for 1 < b < 2, (26) 

Q^-G) for b >2, 

where Q(x) is the unit step function, equal to unity for x > 0 and zero otherwise. The above 
predictions for the tails may be checked against numerically computed Prob(G) using Eq. 
Cl. We show in Fig. [5] a comparison between numerical results and our predictions for 
three representative values of b for a = 0.375. 

Invoking the results on stable distributions discussed above, we may deduce the behavior 
of Prob(p*) for different range of values of b by using Eqs. (TT71) and ([26]) . 

(i) 0 < b < 1: Here, Prob(p*) is a Levy-stable distribution with a tail decaying as a power 
law with exponent (1 + a) for values of a in the range 0 < a < 2 and is a Gaussian 
distribution for a > 2. The mean is finite for a > 1. 

(ii) 1 < b < 2: In this case, Prob(p*) is a Levy-stable distribution with a tail decaying as a 
power law with exponent (1 + v) for values of v in the range 0 < v < 2, that is, provided 


0.50001 


14 










Prob(p*)/c(L) ^ Prob(p*)/c(L) 


(a) 


b=4.0,a=0.75 


b=1.5,a=0.75 (b) 



Figure 6. Scaling plots for the probability distribution Prob(p*) for four representative 
values of the set (b,a), with p* computed using Eq. 1171) . Here, the scaling factors are 
c(L) = Prob max (p*), the maximum value of the distribution, and a(L ) the corresponding 
value of p*. (a) and (b) stand for cases where the system supports condensate formation 
(see condition (ED), and p* coincides with the critical density p c . (c) and (d) denote 
cases where the system does not support condensation. The black lines denote analytical 
predictions, namely, (i) for b = 1.5, a = 0.75 (giving v = 1.5), the inverse Fourier 
transform of the function exp[—|cofc| !/ (l — i/3 sgn(fc)$)], with $ = tan(7Ti//2), /? = 0.995, 
Co = 0.28, which has been shifted by an amount equal to 0.39 to fit the data, (ii) for 
b = 4.0, a = 0.75 and b = 1.5, a = 6.0, a Gaussian with zero mean and variance equal 
to l/(27r), and (iii) for b = 1.5, a = 0.25 (giving u = 0.5), the so-called Levy distribution 
f(x ) = ^ co/(27t) exp[— co/(2(a; — p))\/{x — /r) 3 / 2 ; cq = 0.45 and p = —0.175. 
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a < 4 — 2b. On the other hand, for a > 4 — 2b, the distribution is a Gaussian. The 
mean p* is finite for b + a > 2. 

(iii) b > 2: In this regime, Prob(G) has a finite variance, implying that Prob(p*) is a 
Gaussian; the mean is of course finite. 


Thus, the condition b + a > 2 ensures a finite value of the mean p*, a condition we derived 
earlier, see Eq. (j23|) . based on an analysis of the disorder-average (p*) defined in Eq. (fT8j) . 

As a function of the system size L , the probability distribution Prob(p*) has the scaling 
form (see Fig. |6|): 


Prob(p*) ~ c(L)Q 



o(L))c(L)], 


(27) 


where c(L) is the maximum value of the distribution, while a(L) is the corresponding value 
of p*. The scaling function G(x) has either (a) a Gaussian form, in which case a(L ) is 
independent of L, and c(L) Vl, or, (b) a Levy-stable form, in which case we have 


(i) for 0 < b < 1: a(L) ~ L 1//a 1 and c(L) ~ L 1 1/,Q for 0 < a < 1, and a(L) independent 
of L and c(L ) ~ L 1_1 / a for a > 2, 

(ii) for 1 < b < 2: a(L ) ~ L l / V ~ l and c(L ) ~ L 1_1 G for 0 < v < 1, and a(L ) independent of 
L and c(L) ~ A 1 ” 1 / 1 " for v > 2. 


4. Numerical studies 

In this section, we check our predictions on the existence of a condensate in the RC-ZRP 
by reporting on results obtained by a direct numerical sampling of its canonical steady state 
measure for a given realization of the disorder. When relevant, we also perform Monte Carlo 
(MC) simulations of the RC-ZRP dynamics while starting from the steady state. 

The steady state is generated according to the following algorithm. For a given system 
size L and disorder realization {/*}, a configuration of the system corresponding to a total N 
particles is generated by occupying L — 1 sites independently with n l (0 < n, < /,) particles 
with respective weights /j(fo|/i); here, i = 1,2,... ,L — 1. The deficit number of particles, 
nd = N — Y^i=i 1 n ii w hen positive, is accommodated on the last remaining site % = L with 
the weight If the deficit is negative, the configuration is rejected, and the process 

is repeated all over again. A configuration so generated is run for a typical “equilibration” 
time of order L before performing any analysis of the data. 
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Figure 7. For a fixed disorder realization, (a) shows the occupancy distribution for the full 
system in the steady state, with densities smaller than, equal to, and larger than the critical 
density p c \ the density values are 1.0, 2.28 and 4.0, respectively. The system size is L = 256, 
and we have taken a = 0.5, & = 1.75. The data are obtained by sampling the canonical 
steady state measure using the algorithm detailed in the text. Below p c , the distribution 
shows an exponential decay at large n, while at p c , one has a power-law decay at large n 
with exponent (b + a ), see Eqs. (fl2l) and (l25l) . For p > p c , in addition to the power-law 
behavior, a bump indicating the presence of a condensate appears. The black line stands 
for the power-law behavior ~ n -(b+«)_ (b) shows for p > p c the occupancy distribution in 
the steady state for the full system as well as for sites with the largest, the second largest 
and the third largest capacity. Besides a bump at large n that implies the presence of a 
condensate, one has a power-law decay at small n , with exponent (b+a) for the full system, 
and with exponent b for individual sites, see Eq. (1241) . 
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b=1.25,a=0.5,L=256 



Figure 8. For the same disorder realization, system size and a as in Fig. [7] but for a value 
of b that does not satisfy the conditions (ll6l) to observe condensation, the figure shows the 
occupancy distribution for the full system in the steady state at several densities. The data 
are obtained by sampling the canonical steady state measure using the algorithm detailed 
in the text. In contrast to Fig. [T^a), the distribution at all densities has an exponential 
decay at large n, and there is no extra peak appearing at high densities that corresponds to 
a condensate. Here, we have taken b = 1.25. 


Following the above procedure for parameter values b = 1.75, ct = 0.5 that satisfy the 
conditions (ITU]) to observe condensation, Fig. 17(a) shows the results for the occupancy 
distribution p(n) for the full system with densities below, at, and above the corresponding 
critical density p c , which is computed numerically from Eq. (ITT]) . Consistent with the 
predictions of Section [3J, we find that a distribution that decays exponentially for p < p c 
goes over to one decaying as a power law at p — p c , which at higher densities develops 
an additional bump corresponding to the formation of a condensate. The power-law decay 
exponent equals {b + a), as predicted in Eq. (125]) . Figure [7(a) is to be contrasted with 
Fig. [8] obtained for the same disorder realization, system size and a, but for a value of b 
that does not satisfy the conditions (1T6]) to observe condensation; the distribution decays 
exponentially at all densities. 
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b=1.75,a=0.5,L=256 
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Figure 9. For the same disorder realization and other parameters as in Fig. [71 the figure 
shows the occupancy at sites with the largest, the second largest and the third largest 
capacity, all plotted together as a function of time for a single dynamical evolution of the 
system while starting from the steady state. The data are obtained by performing Monte 
Carlo simulations of the dynamics starting from a steady state configuration. 

4-1. Condensate relocation 

In Fig. [TJ^b) , we contrast the single-site occupancy distribution Pi(n\li) for large-capacity 
sites with the occupancy distribution p{n) for the full system at a density p > p c . The 
single-site distribution is observed to be the same for the site with the largest capacity (let 
us denote it by R) and for the ones with the second and the third largest capacity (denoted 
respectively by indices z 2 and i 3 ). To understand such a behavior, we show in Fig. [9] the 
results of a MC simulation of the dynamics for the same set of parameter values and the 
same disorder realization as in Fig. [0 The occupancy at sites ii, ^3 have been plotted 
as a function of time for a single dynamical evolution of the system starting from a steady 
state configuration. It is evident from the figure that a dip in n± 3 from a value of 0(L ) to 
a value of order 1 is followed by a rise within a short time in either Uj 2 or n l3 from a value 
of order 1 to a value of order L. This implies that the condensate occupies a single site at 
almost all times, but does move between certain sites with a relatively small relocation time. 
The fact that for a given L and a given disorder realization, there are only a finite number 
A r(L,{k}) L 1 “ of sites that have capacities equal to or larger than L implies that the 
condensate can relocate only on this finite subset of sites. A similar relocation dynamics 
of the condensate on a set of sites whose size grows subextensively with L was observed in 
a disordered version of the ZRP studied in [29], in which the disorder enters through hop 
rates. Such a relocation of the condensate on a subset of sites of subextensive size may be 
contrasted with the situation in the homogeneous ZRP where the condensate can relocate 


19 































on any of the other L — 1 sites [361137] . In our case, when the condensate has relocated away 
from one of the sites of this subset to another, the occupancy and fluctuations on the first 
site become identical to ones in the background that did not contain the condensate. This 
explains why the single-site distribution for sites 21 , 22,23 are the same. A detailed analysis 
of the condensate relocation dynamics will be published elsewhere. 

5. Conclusions 

In this paper, we studied a quenched disordered version of the zero-range process (ZRP), a 
nonequilibrium system of particles undergoing biased hopping on a one-dimensional periodic 
lattice. In the model studied, which we refer to as the random capacity zero-range process 
(RC-ZRP), each site has a finite capacity whereby it can hold only a finite number of particles; 
we chose the capacities randomly from the Pareto distribution. We obtained the conditions 
for condensate formation in the RC-ZRP, which derive from an interplay of the capacity 
distribution with the hop rate. In terms of the power-law exponents a > 0 and b > 0 
that characterize respectively the capacity distribution and the hop rate, we derived explicit 
conditions for condensation, namely, b + a > 2 and a < 1. Further, we addressed the sarnple- 
to-sample variation of the critical density to observe condensation, and demonstrated that 
the corresponding distribution is either a Gaussian or a Levy-stable distribution. 

Let us remark on the possibility of observing condensation in the RC-ZRP for 
generalizations of the hop rate 22 ( 22 ) that we studied. Consider, e.g., the choice 22 ( 22 ) = 
1 + b/n a , with a > 0. For a < 1, the function G{1) in Eq. (TT71) converges asymptotically 
to a finite constant for all values of b, yielding a finite p*. As a result, the system supports 
condensation for all values of b, provided a < 1, a condition that derives from the desired 
scaling of l mSLX with system size L. On the other hand, for a > 1, the function G(l ) 
diverges asymptotically for all values of b, so that p * is infinite, and consequently, there 
is no condensate formation in the system. 

We sign off by mentioning a possible follow-up of this work. It would be of interest to 
study the RC-ZRP dynamics in the steady state, and investigate the behavior of time- 
dependent correlation functions. In this regard, a pertinent issue is to address if and 
how quenched disorder manifests itself in the behavior of the dynamic universality class at 
criticality and in the dynamics of condensate relocation, both of which may show significant 
differences from the homogeneous ZRP [[38] . 
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7. Appendix: Characterizing the site capacities — The sum and the maximum 

In this appendix, we summarize some features associated with the capacity distribution (jT]) 
that are relevant to the understanding of the condensation phenomenon in the RC-ZRP 
discussed in the main text. 

Let us start with discussing the behavior of the mean and the variance of the distribution: 
they are both finite for a > 2 and both infinite for 0 < a < 1. In the intermediate 
regime 1 < a < 2, the mean is finite while the variance is infinite. For values of a in 
the range 0 < a < 2, the distribution © is Levy-stable: a linear combination of two 
independently sampled values of / has a distribution identical to P(Z), up to location and 
scale parameters |33ll3o) . For 0 < a < 2, a Levy-stable distribution is characterized by a 
power-law tail with exponent — (1 + a); the distribution is a Gaussian for a = 2. 



200 1000 10000 

L 


Figure 10. Illustrating the validity of the scaling behavior (IdOl) of the disorder-averaged 
ratio (tV max /i max ) for two values of a , one smaller and one larger than 1. The data for 
a = 1.5 have been scaled down by a factor of 4 for convenience of display. 


Now, let us discuss the scaling with system size L of the largest capacity Z max and the 
largest possible number of particles iV max that can be accommodated in the system. For a 
given L and a given realization {U} of the disorder, we have Z max = Max[/i,/ 2 , • • •, Il\ and 
iV max = h- Since the Z*’s are sampled independently from the common distribution (JTJ), 
the probability distribution of / max is 

f Umax . L—1 

Prob(Z max ) = LP(Z max ) (^ dl P(l)J . (28) 
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In the limit L —> oo, the distribution Prob(/ max ) decays for large l max S> 1 as Prob(/ max ) 
exp(—L/“" x ), which implies the following scaling of / max with L: 

/max ~ L 1/a , (29) 


valid for all values of a > 0. 

As to the behavior of 7V max , for a > 1 , when P{1) has a finite mean, one may apply the 
law of large numbers to deduce that N mgx = aL in the limit L —* oo, with a = f dl IP{1 ) 
finite. For 0 < a < 1, on the other hand, the mean of P{1) is infinite, the law of large 
numbers breaks down, and A r max is dominated by contributions from capacities of order Z max . 
Thus, we anticipate N max ~ L 1 11 /", which would imply 


N max \ I 1 for 0 < a < 1, 

lmax ! | L l ~ 1,a for a > 1. 


(30) 


Figure [10] illustrates the validity of the above scaling for representative values of a. In fact, 
the full distribution of the ratio N max /l max is known (see (31], page 465), which leads to 


N n 


a 


, 0 < a < 1. (31) 

^max ' ^ 

The above result is confirmed in Fig. [TO] Note that Eq. (130|) suggests that there are several 
sites other than the site with capacity l max which have capacities of order L l ^ a . 
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